Non-Markov dynamics and phonon decoherence of a double quantum dot charge qubit 
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Abstract 

In this paper we investigate decoherence times of a double quantum dot (DQD) charge qubit due to it coupling with acoustic 
phonon baths. We individually consider the acoustic piezoelectric as well as deformation coupling phonon baths in the qubit 
environment. The decoherence times are calculated with two kinds of methods. One of them is based on the qusiadiabatic 
propagator path integral (QUAPI) and the other is based on Bloch equations, and two kinds of results are compared. It is 
shown that the theoretical decoherence times of the DQD charge qubit are shorter than the experimental reported results. It 
implies that the phonon couplings to the qubit play a subordinate role, resulting in the decoherence of the qubit. 

PACS numbers: 73.63.Kv, 03.65.Yz, 03.67.Lx 



I. INTRODUCTION 



Solid state qubits are considered to be promising can- 
didates for realizing building blocks of quantum informa- 
tion processors because they can be scaled up to large 
numbers. The double quantum dot (DQD) [J H Hll U 
charge qubit is one of these qubits. Two low-energy 
charge states are used as the local states |0) and |1) in the 
qubit. The qubit can be controlled directly via external 
voltage sources. There are some effective schemes to pre- 
pare the initial states and read out the final states of the 
qubit [6( . So it is considered that decoherence may be the 
central impediment for the qubit to be taken as the cell 
of quantum computer. Finding out the primary origin or 
the dominating mechanism of decoherence for the qubit 
is a basal task for overcoming the impediment. Experi- 
mentally, many attempts 0, for detecting the decoher- 
ence time of this kind of qubit have been performed. The 
decoherence has also been investigated theoretically. In 
2000, Fedichkin et al. @ investigated the Born-Markov 
type electron-phonon decoherence at large times due to 
spontaneous phonon emission of the quantum dot charge 
qubits. Recently, Vorojtsov et al. 0] studied the deco- 
herence of the DQD charge qubit by using Born-Markov 
approximation. But, as it has been pointed out, the use 
of the Born-Markov approximation is inappropriate at 
large tunneling amplitudes. The method is expected to 
become increasingly unreliable at DQD with larger inter- 
dot tunneling amplitudes. Wu et al. [ll| investigated the 
decoherence in terms of a perturbation treatment based 
on a unitary transformation. The Born-Markov approx- 
imation has not been used in the method but it neglects 
some terms of the effective Hamiltonian with high excited 
states. This kind of processing introduces an new approx- 
imation which has not been estimated to the affects of 
the dynamics. Fedichkin et al. [T2I IT^ | studied the error 
rate of DQD charge qubit with short-time approximation. 
This method is accurate enough in adequate short time. 
But the decoherence in a moderately long time is also in- 
teresting. Recently, Thorwart et al. investigated the 
decoherence of the DQD charge qubit in a longer time 
with a numerically exact iterative quasiadiabatic propa- 



gator path integral (QUAPI) This method is proved 
valid in investigating the qubit decoherence jig . In Ref. 
[lH , Thorwart et al. considered the coupling of longitu- 
dinal piezoelectic acoustic phonons with the investigated 
qubit and neglected the contribution of the deformation 
acoustic phonons to decoherence. These two kinds of 
phonons may constitute two kinds of different coupling 
baths in the environment of the qubit. We call the former 
the piezoelectric coupling phonon bath (PCPB) and the 
latter the deformation coupling phonon bath (DCPB). 
Comparing Thorwart 's result and the reported experi- 
mental value they found that the theory predicts the de- 
coherence time of the DQD charge qubit is two orders 
of magnitudes smaller than the experimental one. Thus, 
Thorwart et al. conclude that the piezoelectric coupling 
phonon decoherence is a subordinate mechanism in de- 
coherence of the DQD charge qubit. Recently, Wu et 
al. gave the spectral density functions of PCPB as 
well as DCPB. Then how about the DCPB to the de- 
coherence of the DQD charge qubit? In other words, 
is the decoherence of the DQD charge qubit induced by 
DCPB also subordinate? In this paper we shall use an 
iterative tensor multiplication (ITM) ^jjj scheme derived 
from the QUAPI to study the decoherence times of the 
DQD charge qubit not only in PCPB but also in DCPB. 
In order to validate if our result is in accordance with 
Thorwart 's result we at first investigate the decoherence 
times of the qubit in PCPB. Then we shall investigate the 
decoherence times of the qubit in another bath, DCPB, 
which will show that the influence of the DCPB to the 
decoherence of the DQD charge qubit is also subordinate 
because it results in a shorter decoherence time than the 
experimental value of 1 ns 0, ■ 



II. MODELS 

The DQD charge qubit consists of left and right dots 
connected through an intcrdot tunneling barrier. Due to 
Coulomb blockade, at most one excess electron is allowed 
to occupy the left and right dot, which defines two basis 
vectors |0) and |1) . The energy difference e between these 
two states can be controlled by the source-drain voltage. 
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Neglecting the higher order tunneling between leads and 
the dots, the effective Hamiltonian in the manipulation 
process reads 0, Q] 



{M q b\ 



H eff = hT c a x + h}^ojgblb q + ha z } j i 

Q 1 

(1) 

Here, T c is the interdot tunneling, a x and <r z are Pauli 
matrix, b q (b q ) are the creation (annihilation) opera- 
tors of phonons, tkj q is the energy of the phonons, and 
M q = Cqj yj2m q uj q h where C q are the classical coupling 
constants of qubit-phonon system. We call the collec- 
tive coupling phonons to the qubit in the environment 
a phonon bath. In order to obtain the reduced density 
matrix of the qubit in the system, one should know the 
coupling coefficients M q , but in fact we need not know 
the details of each M q because all characteristics of the 
bath pertaining to the dynamics of the observable system 
are captured in the spectral density function |l7l Il8| 



J(«) =J2\M q \ 2 8(u>-u> q ). 



(2) 



It is pointed out that the spectral density of PCPB is 

JP* (uj) = g oj (\-^± sin —\ e~^? . (3) 
\ oj UJ d J 

Here, u>d = s/d and lui = s/l, where d denotes the center- 
to-center distance of two dots, I the dot size, s the sound 
velocity in the crystal, and 
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Here, M is the piezoconstant, g is the density of the 
crystal, and x is the rate of transverse to the longitudinal 
of sound velocity in the cr yst al, (see for example Refs. 

[HE!). As in Refs - dill in this P a P er we set the 
sound velocity in the GaAs crystal s ~ 5 x 10 3 m/ 



n/s . 



With the parameters of GaAs in Ref. 19], Wu et 
propose a value g pz w 0.035 (ps)~ 2 . The spectral density 
of DCPB is 



9df0 



where 



— sin — 
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Here, S is the deformation potential. In the same paper, 
Wu et al. also propose a value gdf ~ 0.029 (ps)~ 2 . One 
can investigate the dynamics and then the decoherence 
of the open qubit with the help of the definite spectral 
density functions of the baths. Before investigations of 
decoherence of the DQD charge qubit we introduce an op- 
timal numerical path integral method, the ITM method 
in the following section. 



III. QUAPI AND ITM 

In the following, we firstly review the QUAPI and then 
the ITM scheme. Suppose the initial state of the 
qubit-bath system has the form 



R(P)=p(0)® Pbath (0), 



(5) 



where p (0) and p oat h (0) are the initial states of the qubit 
and bath. The evolution of the reduced density operator 
of the open qubit 

p (s", a'; t) = Tr bath (s"\ e^ Ht / h p (0)® Pbath (0) e im ' h \s') , 

(6) 

is given by 

p(s",s';t) 

= E E ••• E E E ••• E 

4=±ls+=±l s+_ 1 =±lso=±ls-=±l s-_ 1= ±l 
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X <4|p(°) | s o) 
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where the influence functional is 



Tr, 



bath 



^-iH snl ,(s")At/2h -i-ff c „„(s+ _ 1 )At/2H 



_ , . e -iH env (s+)At/2H /qN e iH anv (s Q )At/2h 

JHenv (s K _j) At /2S iff e „„ (s/)At/2R 



(8) 



Here, Hq is a reference Hamiltonian that in general 
depends on the coordinate and momentum of the 
system. In the qubit system, it usually depends on 
Pauli matrixes a x and a z . The H env is defined as 
H env = H — Hq. In our system we set Hq — hT c a x . 
The discrete path integral representation of the qubit 
density matrix contains temporal nonlocal terms 

1 (4.4. • • ') s jv-i> s "> s o> s i> ■ ■ ■> s iV-i> s '; A< ) which 
denotes the process being non-Markovian. With the 
quasiadiabatic discretization of the path integral, the 
influence functional, Eq.JHJ) takes the form 



/ = exp 



N k 

■EE' 

■ k=0 k'=0 



(9) 

where sj- ~ s" and s N = s' . The coefficients r] kk , 
can be obtained by substituting the discrete path into 
the Feynman- Vernon expression. Their expressions have 
been shown in Ref. flS ]. Thus, the influence functional 
can be expressed with a product of terms corresponding 
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to different Ak as 



N N-l JV-Afe 

i = uio(4) n h(4,4 + i) n w 

k=0 k=0 k=0 



k ' °k+Ak) 



N-Ak 

n 



(10) 



Here, A/c = k — k', where fc' and fc are points of discrete 
path integral expressions, (see Ref. [l^) ar± d 



h (-?) 

^Afc (sf ,S.f +Afe ) 



exp 



exp 



+Afe ^i+Afc; 



X (Vi+Ak,i s i ~ Vi+Ak,i S i )} , Afc > 1. 

(11) 



The length of the memory of the time can be estimated 
by the following bath response function 



it) 
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dcoJ x (uj) 



coth cos tot — % sin 

2 



(12) 

Here, the superscript x denotes the bath type, [3 — 
X/ksT, where ks is the Boltzmann constant, and T is 
the temperature. It is shown that when the real and 
imaginary parts behave as the delta function 5 (t) and 
its derivative 6 (t) , the dynamics of the reduced density 
matrix is Markovian. However, if the real and imaginary 
parts are broader than the delta function, the dynam- 
ics is non-Markovian. The broader the Relet 21 (t)} and 
Imla^ (t)] are, the longer of the memory time will be. The 
broader the Relet 21 (t)} and Imfa 21 (t)] are, the more seri- 
ous the Markov approximation will distort the practical 
dynamics. In Fig.l we plot the Re[a pz (t)] and Im[ct pz (t)] 
of the PCPB and in Fig.2 we plot the Re[a rf/ (*)] and 
Im[a rf/ (*)] of the DCPB. 
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rapidly as the "interaction distance" increases. In the 
ITM scheme the interaction can be taken into account at 
each iteration step. The reduced density matrix at time 
t = NAt (N even) is given as 



~p(s%,NAt)=A^ (s 



where 



A« (*± +1 ;(* + l)Ai) 



%;NAt)l (s±) 



iA® (sf;kAt) . 



Here, 



T (2Afe max ) / ± ± ± 1 

fc+Afc max -l 

n K(si,si +1 )i {st)h(st,st +1 ) 

n—k 

x/ 2 {st,sf l+2 ) ...I A /c max («n' s n+AAw) ' 

and 
where 



-i,o) = (4\p. (°)k) 



K i4>4i-i) = {4 + i\ c M-iH At/h)\s+) 
x (s^\enp(iH At/h) |s^ +1 ) . 

In the ITM scheme a short-time approximation instead 
of the Markov approximation is used. The approxima- 
tion makes an error of the short-time propagator in order 
(At) 3 , which is small enough as we set the time step At 
very small. It is shown that when the time step At is 
not larger than the characteristic time of the qubit sys- 
tem, which can be calculated with 1/T C , the calculation 
is accurate enough [2(| • In particular, the scheme does 
not discard the memory of the temporal evolution, which 
may be appropriate to solve the decoherence of qubit. In 
the following section we shall use the ITM scheme to 
study the decoherence times of the DQD charge qubit in 
PCPB and DCPB. 



Fig.l 



We see that the memory times are about r^ em = 
1 x 1CT 11 s for PCPB and r# em = 2 x KT 11 s for 
DCPB {where the points beyond ±1 x 10 -11 s have not 
been plotted for clearly distinguishing the Re[a * (t)] and 
Imfa* (t)] in the same figure}. Due to the nonlocality, 
it is impossible to calculate the reduced density matrix 
by Eq.Q in the matrix multiplication scheme. However, 
the short range nonlocality of the influence functional 
implies that the effects of the nonlocality should drop off 



IV. DECOHERENCE OF DQD CHARGE QUBIT 

To measure effects of decoherence one can use the en- 
tropy, the first entropy, and many other measures, such 
as maximal deviation norm, etc. (see for example Ref. 
ppf). However, essentially, the decoherence of a open 
quantum system is reflected through the decays of the 
off-diagonal coherent terms of its reduced density matrix. 
The decoherence is in general produced due to the inter- 
action of the quantum system with other systems with a 
large number of degrees of freedom, for example the de- 
vices of the measurement or environment. The decoher- 
ence time denoted by measures the time of the initial 
coherent terms to their 1/e times, namely, p i {n,m) 
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Pf (n, m) = pn (n, m) /e. Here, n ^ m, and n, m = or 
1 for qubits. In this paper, we investigate the decoher- 
ence times via directly describing the evolutions of the 
off-diagonal coherent terms, instead of using any mea- 
sure of decoherence. In our following investigations, we 
suppose the temperature T = 30 mK and the cut-off fre- 
quency of the bath modes uic = 5 (ps) -1 . We set the ini- 
tial state of the qubit to p(0) = | (|0) + |1)) ((0| + (1|) , 
which is a pure state and it has the maximum coher- 
ent terms, and the initial state of the environment is 
Pbath (0) = Uk e ~ 0Mk / TT k (e-P M *) . In the calculations 
we set LOd = 0.02 (ps) -1 , T c — O.Iuji according to Ref. 
[Tl|. and two kinds of cases = 0.5 (ps) -1 and = 0.7 
(ps) -1 are calculated. 

Decoherence time obtained from ITM scheme: In the 
following, at first, we use the ITM scheme investigating 
the decoherence time of the DQD charge qubit. The 
evolutions of the coherent elements of the reduced density 
matrix of the DQD charge qubit in PCPB and DCPB 
are plotted in Figs. 3 and 4. Here, we simply choose 
Afc max = 1 and At = 1 x 10 -11 s for PCPB and At = 2 x 
10 -11 s for DCPB in the ITM scheme. These choices of 
the time steps are feasible as we consider that it should be 
not smaller than the memory times of the baths, because 
the latter is about T^ em w 1 x 10 -11 s for PCPB and 
r^ em « 2 x 10" 11 s for DCPB (see Figs. 1 and 2). It is 
also appropriate as we consider that the time steps should 
not be larger than the characteristic time of the qubit, 
because the characteristic time of the qubit is about 2 x 
10- n s. 



Fig.3, 



FigA 



Helped with detailed numerical analyses, we can obtain 
that the decoherence times of the DQD charge qubit in 
PCPB are about r^ 2 ~ 97 ps [when u>i = 0.7 (ps) -1 ] and 
t V 2 « 118 ps [when uji = 0.5 (ps) -1 ]. Similarly, we can 
obtain that the decoherence times of this qubit in DCPB 
are about t 2 = 1-04 ps [when u)\ = 0.7 (ps) -1 ] and 
t 2 = 3.5 ps [when usi = 0.5 (ps) -1 ]. It is shown that the 
DCPB behaves more destructively than the PCPB does 
to the coherence of the DQD charge qubit. A further 
calculation shows that the decoherence time will increase 
with the decreasing of T c . 

Decoherence time calculated on Bloch equations: It is 
well known that the decoherence time can be calculated 
based on Bloch equations. In the following, we calcu- 
late the decoherence time of the DQD charge qubit in 
PCPB and DCPB with the Bloch equation method. In 
this method, the relaxation and dephasing times can be 



evaluated from the spin-bosonic model with Bloch equa- 
tions 0, 0| . For our model, they are [2^] 

t -1 = t -1 = i J K) coth (phuo/2) , 

where uq = 2T C is the natural frequency of the DQD 
charge qubit. By using the parameters of the DQD 
charge qubit and PCPB bath as above, we can calculate 
the decoherence times with this method as t v ^ « 122.3 
ps [when lq{ = 0.7 (ps) 1 ] and t^ z » 192.2 ps [when 
u>l = 0.5 (ps)~ ]. Similarly, we can obtain the decoher- 
ence times of the DQD charge qubit in the DCPB with 
this method as ~ 3.18 ps [when uii = 0.7 (ps) ] and 
t 2 ~ 12.6 ps [when ui = 0.5 (ps) ]. It is shown that 
the decoherence times obtained from the ITM scheme 
are shorter than those obtained on Bloch equations. We 
suggest that the differences are derived from the different 
choices of approximation schemes. The Bloch equations 
are in general derived from the Markov approximation 
which discards the memory of baths in the derivation 
of dynamical evolution. The decoherence of the qubit 
obtained on Bloch equations is similar to the "resonant 
decoherence" [2lJ obtained from the Fermi golden rule. It 
is not accurately equal to the actual decoherence except 
that the "nonresonant decoherence" very small. 

Decoherence time derived from the quality factor: We 
like to compare our results obtained from the ITM 
scheme based on QUAPI with Thorwart's results which 
are also obtained from QUAPI. Thorwart et al. in- 
vestigated the PCPB case and they obtained the qual- 
ity factor instead of the decoherence time. By using 
a set of parameters of the DQD charge qubit and the 
PCPB they obtained the quality factor of the qubit 
as Q pz — 336, which corresponds to decoherence time 
t P 2 = Q pz ir /w' pz ~ 115.9 ps, where u>' = ujq + Aw, and 
Alo is the bath-induced shift I22| in the natural frequency 
uj = 2T C . From Fig.l of Ref. |ll| we see Auj pz w 1.75u c 
and Audf ~ 1.65o; c . Their used parameters [T = 10 mfc, 
T c k 0.07 (ps) ] have a little difference from ours. But 
we have calculated that the difference does not result 
in much decoherence time departure. It is meant that 
our results is in accordance with Thorwart's result. On 
the other hand, from our decoherence time r 2 ~ 3.5 ps 
of the qubit in DCPB we can obtain its quality factor 
Qdf = rfu'df/n w 8. 

V. DISCUSSIONS AND CONCLUSIONS 

In this paper we investigated the decoherence times 
of the DQD charge qubit in PCPB and DCPB with the 
ITM scheme based on QUAPI. The decoherence times 
are also calculated based on Bloch equations. The results 
derived from the two kinds of methods are compared to 
each other. It is shown that the latter are longer than the 
former. We think this results from the different choices 
of approximation schemes because the Markov approxi- 
mation used in the latter method discards the memory of 
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the baths. On the other hand, Hayashi et al. 0,0 have 
detected that the decoherence time of the DQD charge 
qubit is about 1 ns as T c ~ 0.07 (ps)' 1 . The exact ITM 
theoretical decoherence times are two orders of magni- 
tude and five orders of magnitude smaller than the ex- 
perimental value even when we consider the DQD charge 
qubit in independent PCPB and DCPB. These can fi- 
nally and without accident lead to the conclusion that 
the phonon decoherence is a subordinate mechanism in 
the DQD charge qubit. In general, besides the phonon 
couplings' decoherence, the qubit can also result in deco- 
herence from electromagnetic fluctuations (with Ohmic 



noise spectrum), cotunneling effect, background charge 
fluctuations (with 1// noise spectrum), and so on. To 
find out the dominating mechanism of the DQD charge 
qubit decoherence and to the best of our abilities to sup- 
press the central decoherence resources are important 
challenges in the quantum computation field. 
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VI. FIGURE CAPTIONS 

Fig.l: Real (line) and imaginary (short lines) parts 
of the response function of the piezoelectric coupling 
phonon bath (PCPB). Here, we set the temperature 
T = 30 mK, and uj d = 0.02 (ps)" 1 , uji = 0.5 (ps)~\ 
g pz = 0.035 (ps)~ 2 . 

Fig. 2: Real (line) and imaginary (short lines) parts of 
the response function of the deformation coupling phonon 
bath (DCPB). Here, we set g df = 0.029 (ps)~ 2 , and other 
parameters are same as those in Fig.l. 

Fig. 3: The evolutions of the off-diagonal elements of 
the reduced density matrix for the DQD charge qubit in 
PCPB when loi = 0.5 (ps^ 1 (line) and uji = 0.7 (ps)' 1 
(short lines). Here, the cutoff frequency is uj c — 5 (ps)^ 1 , 
other parameters are same as those in Fig.l. The initial 
state of the qubit and environment are described in the 
text. 

Fig. 4: The evolutions of the off-diagonal elements of 
the reduced density matrix for the DQD charge qubit in 
DCPB when cui = 0.5 (ps)- 1 (line) and cj ; = 0.7 (ps)" 1 
(short lines). Here, the parameters are same as those in 
Figs.2 and 3. 



5 



4x10 



3x10 



2x10 



1x10 



■1x10 



-2x10 20 H 



-3x10 



-4x10 




Xian-Ting Liang 
PRB-BG 10449 
Fig.1 



-1.0x10~ 11 -5.0x1 0" 12 



0.0 



5.0x1 0" 12 1.0x10" 11 



t(s) 



1.5x10 

O 1.0x10 
Q 



46 



46 



C 

o 

o 
c 

Z5 

M— 

CD 
C/3 
C 

o 

Q. 

W 
CD 



5.0x10 



45 



0.0 



-5.0x10 



45 



q: -1.0x10 



46 



-1.5x10 



46 



Re[a? z (t)] 
lm[gP z (t)] 




-1.0x10" 11 -5.0x10" 12 0.0 



5.0x1 0" 12 1.0x10" 11 



Xian-Ting Liang 
PRB-BG10449 
Fig.2 



t(s) 



0.5- 

0.4- 

0.3- 

0.2 

0.1 

0.0-I 








co=0.5 (ps) 



-1 



— 03=0.7 (ps) 



10 



20 

t(ps) 



30 



40 



Xian-Ting Liang 
PRB--BG 10449 
Fig.4 



